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CN : ABSTRACT 

. Context. The interaction of plasma motions and magnetic fields is an important mechanism, which drives solar activity in all its 

facets. For example, photospheric flows are responsible for the advection of magnetic flux, the redistribution of flux during the decay 
of sunspots, and the built-up of magnetic shear in flaring active regions. 

Aims. Systematic studies based on G-band data from the Japanese Hinode mission provide the means to gather statistical properties 
of horizontal flow fields. This facilitates comparative studies of solar features, e.g., G-band bright points, magnetic knots, pores, and 
sunspots at various stages of evolution and in distinct magnetic environments, thus, enhancing our understanding of the dynamic Sun. 
Mettiods. We adapted Local Correlation Tracking (LCT) to measure horizontal flow fields based on G-band images obtained with 
, the Solar Optical Telescope on board Hinode. In total about 200 time-series with a duration between 1-16 h and a cadence between 

15-90 s were analyzed. Selecting both a high-cadence (A^ = 15 s) and a long-duration (AT = 16 h) time-series enabled us to optimize 
and validate the LCT input parameters, hence, ensuring a robust, reliable, uniform, and accurate processing of a huge data volume. 
' Results. The LCT algorithm produces best results for G-band images having a cadence of 60-90 s. If the cadence is lower, the velocity 

^1^, of slowly moving features will not be reliably detected. If the cadence is higher, the scene on the Sun will have evolved too much to 

bear any resemblance with the earlier situation. Consequently, in both instances horizontal proper motions are underestimated. The 
most reliable and yet detailed flow maps are produced using a Gaussian kernel with a size of 2560 km x 2560 km and a full-width-at- 
half-maximum (FWHM) of 1200 km (corresponding to the size of a typical granule) as sampling window. 
^ I Conclusions. Horizontal flow maps and graphics for visualizing the properties of photospheric flow fields are typical examples for 

^ value-added data products, which can be extracted from solar databases. The results of this study will be made available within the 

'small projects' section of the German Astrophysical Virtual Observatory (GAVO). 
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(N( . 1 . Introduction subsumed under FT techniques, since artificial tracer particles 



are introduced to follow the footprints of local intensity min- 
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^ Data from space do not suffer the deleterious effects of Earth's ima. The ball tracking method works well for granulation so that 
turbulent atmosphere, which blur and distort images so that fea- it is a good choice for the characterization of supergranulation 
tures might fade into obscuration, making it difficult to follow (Potts & Diver, 2008). Nevertheless, it might introduce a scale 
them from image to image. The huge volume of Hinode G-band dependence, when tracking on other (small-scale) features such 
images with good spatial resolution, cadence, and coverage pro- as bright points, penumbral grains, and umbral dots, 
vide time-series of consistent quality to quantify photospheric ^he second class is based on Local Correlation Tracking 
^ proper motions, which can be used m comparative studies. ^lCT), where displacement vectors are derived by cross- 
Various techniques have been developed in the past decades correlating small regions in consecutive images of a time- series, 
to measure horizontal proper motions on the solar surface. The principles of LCT are laid out in the seminal work of 
Basically, they can be divided in two classes. The first class in- November & Simon (1988). Leaving behind the underlying ve- 
cludes Feature Tracking (FT) methods, which follow the foot- locity assumption of LCT Schuck (2006) developed a technique 
prints of individual features in images of a time-series (see e.g. to track optical flows named the differential affine velocity esti- 
Strous, 1995). Tracking facular points of opposite magnetic po- mators (DAVE). Chae & Sakurai (2008) presented a formulation 
larity in an emerging flux region effectively demonstrated the of the non-linear case and called it accordingly the nonlinear 
potential of FT techniques (Strous et al., 1996). However, the affine velocity estimators (NAVE). These authors also provide a 
image processing steps (segmentation, labeling, and identifica- detailed parameter study of LCT, DAVE, and NAVE based on 
tion) rely on prior knowledge about the object under investiga- images reflecting analytical solutions of the continuity equation 
tion. Therefore, FT methods seem to be better suited for case as well as on magnetogram data from Hinode and MHD simu- 
studies rather than bulk processing of huge data volumes, where lations. Similarly, Welsch et al. (2007) compared these velocity 
the reduction of dimensionality is a desirable feature. The ball- inversion techniques based on MHD simulations, where the ve- 
tracking method developed by Potts et al. (2004) can also be locity field is known. 
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Fig.l. Calibrated G-band images of active regions (a) NOAA 10960 and (b) NOAA 10930 observed on 2007 June 4 and 2006 
December 7, respectively. The FOV is 1 1 x 1 1 T'. The annotation of the axes refers to heliocentric coordinates given in seconds 
of arc. The region in the rectangular white box is used for further analysis zooming in on a quiet Sun region containing numerous 
G-band bright points outlining supergranular boundaries. The images were normalized such that the mode of the quiet Sun intensity 
distribution corresponds to unity. 



Various other technical aspects have to be considered while 
implementing techniques for optical flow tracking. Since inter- 
polations are required in the various data reduction steps, we 
refer to Potts et al. (2003), who describe the systematic errors, 
which can be introduced by unsuitable interpolation schemes. 
More recently, Lofdahl (2010) discussed image-shift measure- 
ments in the context of solar wavefront sensors, which are also 
applicable to the LCT implementations. 

Hinode G-band images oflfer an excellent opportunity for 
systematic statistical studies of flow fields, because of their uni- 
form quality in the absence of seeing distortions, thus alowing 
us to directly compare flow fields in diflferent solar settings. In 
Sect. 2, we introduce the high-cadence and long-duration data 
sets, which are used to find the optimal parameters for comput- 
ing LCT maps based on time-series of Hinode G-band images. 
The implementation of the algorithm is described in Sect. 3. 
Section 4 presents the results of this parameter study justify- 
ing our choice of LCT parameters, which will be used to create 
value-added data to complement the Hinode database. The most 
important parameters are summarized in Sect. 5, which also in- 
troduces some of the future uses of such a database. 



2. Observations 

Images obtained in the Fraunhofer G-band (bandhead of the CH 
molecule at /1430.5 nm) have high contrasts, and small-scale 
magnetic features can be simply identified with bright points 
(Berger et al., 1995). Despite the observational advantages of 
such "proxy-magnetometry" (Leenaarts et al., 2006), the theo- 
retical description of the molecular line formation process is far 
from easy (cf. Sanchez Almeida et al., 2001 ; Steiner et al., 2001 ; 
Schiissler et al., 2003). LCT techniques, however, can take full 
advantage of the high contrast and the rich structural contents of 
G-band images. On Hinode (Kosugi et al., 2007) such observa- 
tions are carried out by the Broad-band Filter Imager (BFI) of 
the Solar Optical Telescope (SOT, Tsuneta et al., 2008). 

Our initial selection criteria were that at least 100 G-band 
images had to be recorded on a given day, which should have 
additionally a cadence of better than 100 s. It turned out that 



these criteria restricted us to data with only half the spatial reso- 
lution (O.ir' pixel"^), where 2x2 pixels were binned into one. 
In total 48 data sets with 2048 x 1024 pixels and 153 data sets 
with 1024 X 1024 pixels were selected for further analysis. The 
time intervals covered by these data sets range from one to 16 
hours. The bulk statistical analysis will be presented in forth- 
coming publications. Here, we will discuss in detail the LCT 
algorithm and justify our choice of input parameters. For this 
purpose, we selected two data sets: one with a high cadence and 
another one with a long duration. In addition, we picked a data 
set without binning to study the dependence of flow maps on the 
spatial resolution of the input data. 

2.1. High-cadence sequence 

LCT depends on several input parameters such as the time in- 
terval between successive images and the sampling window's 
size and form. We analyzed an one-hour time-series with a ca- 
dence of 15 s to validate the intrinsic accuracy of the LCT al- 
gorithm. The data were captured from 14:27-15:27 UT on 2007 
June 4 (see Fig. la). The time-series contains 238 images with 
1024 X 1024 pixels. The observations were centered on active 
region NOAA 10960, which was located on the solar disk at he- 
liocentric coordinates E630'' and S125'' (ji = 0.75). The active 
region was in the maximum growth phase and had a complex 
magnetic field configuration. NOAA 10960 was classified as a 
j3y6-YQgion and was the source of many M-class flares including 
a major M8.9 flare at 05:06 UT on 2007 June 4, which has been 
analyzed in a multi- wavelength study by Kumar et al. (2010). 

2.2. Long-duration sequence 

Solar features evolve on diflferent time scales from about five 
minutes for granulation to several tens of hours for supergranu- 
lation. Obviously, the time over which individual LCT maps are 
averaged plays a decisive role in the interpretation of such aver- 
age flow maps. Hence, we selected a time-series with 16 hours of 
continuous data captured on 2006 December 7 (see Fig. lb). This 
long-duration sequence starts at 02:30 UT and ends at 18:30 UT. 
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It includes 960 images with 1024 x 1024 pixels and has a ca- 
dence of 60 s. Due to memory constrains imposed by the sub- 
sonic filtering we choose an area of 210 x 595 pixels centered on 
a region with granulation and G-band bright points (white box 
Fig. lb). This region is to the west of active region NOAA 10930 
located at heliocentric coordinates E777''and {ji = 0.59). 
The sunspot group was classified as a jSyS-region exhibiting a 
complex magnetic field topology and producing numerous C-, 
M-, and X-class flares. This region has been extensively studied, 
especially around the time of an X3.4 flare on 2006 December 13 
(e.g., Schrijver et al., 2008). LCT techniques were used by Tan 
et al. (2009) to study horizontal proper motions associated with 
penumbral filaments in a rapidly rotating ^-spot. 

2.3. High-spatial resolution sequence 

Only a few data sets of G-band images exist with the full spa- 
tial resolution of 0.055'' pixel" ^ and a cadence suitable for LCT. 
We selected a one-hour time- series, which was acquired start- 
ing at 04:00 UT on 2006 November 26. This high-spatial res- 
olution sequence contains 118 images with 2048 x 2048 pixels 
and has a time cadence of about 30 s. The observations were 
centered on a quiet Sun region near disk center at heliocentric 
coordinates El 08'' and SI 25" (ji = 0.99) containing only a few 
G-band bright points and no major magnetic flux concentrations. 

3. Implementation of the LCT algorithm 

3.1. Preprocessing of the G-band images 

The data analysis is carried out in the Interactive Data Language 
(IDL)^ Data sets are split in 60-minute sequences with 30 min 
overlap between consecutive sequences. In preparation for the 
LCT algorithm basic data calibration is performed, which con- 
sists of dark current subtraction, correction of gain, and removal 
of spikes caused by high energy particles. Figure 1 contains sam- 
ple G-band images for both the high-cadence and long-duration 
data sets after basic calibration. After initial data calibration, the 
geometric foreshortening is corrected and images are resampled 
in a regular grid with spacing of 80 km x 80 km, i.e., the images 
appear as if observed at the center of the solar disk. Residual 
eflTects of projecting the surface of a sphere onto a plane are ne- 
glected, since the FOV of the G-band images is still relatively 
small. A grid size of 80 km was chosen so that the fine structure 
contents of the G-band images are not diminished. Pixels close 
to the solar limb are projected onto several pixels in planar co- 
ordinates. Thus, the accuracy of the flow maps at these locations 
is not as good as for locations close to disk center. 

In a 60-minute sequence are / = 0, 1, 2, . . . , L - 1 cali- 
brated and deprojected images, where L is the total number 
of images in a particular sequence. For an image with N x M 
pixels the intensity distribution is represented by i(x,y) with 
X = xo, xi, . . . , xn-i and y = yo,yi,- • . , Jm-i as pixel coordi- 
nates. The indices are typically dropped to ease the notation. 
The data processing makes extensive use of the Fast Fourier 
Transform (FFT). The FFT of the intensity distribution i(x,y) 
is simply denoted by 9^(i(x, y)). 

3.2. Aligning the images within a time-series 

In principle, images could be aligned using the pointing infor- 
mation of the spacecraft. However, we calculate shifts between 

^ www.ittvis.com 



consecutive images //-i(x, y) and //(x, y) by computing the cross- 
correlation using only the central part of the images, which is 
half of the original image size. These shifts are then applied in 
succession to align all images with respect to the first image us- 
ing cubic spline interpolation with subpixel accuracy. The signa- 
ture of the 5-minute oscillation is removed from the time-series 
by applying a 3D Fourier filter. This filter, sometimes called a 
subsonic filter, has a cut-oflT velocity of 8 km s"^ corre- 
sponding to the photospheric sound speed. Since the subsonic 
filter uses a 3D Fourier transform some edge eflTects are some- 
times noted for the first and last few images of a time- series. 
For this reason we decided to discard the images during the first 
and last two minutes of the time- series after applying the sub- 
sonic filter. Therefore, the final time-series is shortened by this 
amount of time (see Sect. 3.5). 

3.3. LCT algorithm 

The LCT algorithm is based on ideas put forward by November 
& Simon (1988). The algorithm has been adapted to subim- 
ages with sizes of 32 x 32 pixels corresponding to 2560 km 
X 2560 km, so that structures with dimensions smaller than a 
granule will contribute to the correlation signal. Since cross- 
correlation techniques are sensitive to strong intensity gradients, 
a high-pass filter was applied to the entire image suppressing 
gradients related to structure larger than granules. The high-pass 
filter is implemented as a Gaussian with a FWHM of 15 pixels 
(1200 km). To indicate that we refer to an subimage with 32 x 
32 pixels and not the entire image, we use the notation i(x\y'). 
The Gaussian kernel used in the high-pass filter then becomes 

where cr = FWHM/(2 V21n2) and r(x\y) = {x'^+y^Y'^. The 
high-pass filtered image can be expressed as 

/high(-^, y) = i{x, y) - i(x, y) g(x\y) , (2) 

where denotes a convolution. The result is an image rich in 
detail, where the low spatial frequencies have been removed. 

The core of the LCT algorithm is the cross-correlation 
ci(x,y, x\y') computed over a 32 x 32 pixel region centered on 
the coordinates (x,y) for each pixel in image pairs ii-i(x,y) and 
ii(x, y), which can be written as 

ci(x,y,x\y) = %{r-'[r{ii-i(x,y,x\y)g(x\y))' 

r\ii(x,y,x\y)g(x\y))]]d(x\y) , (3) 

where g(x\y') denotes a weighting function also serving as 
an apodising window. This function has the same form as the 
Gaussian kernel previously used in the high-pass filter. This en- 
sures that the displacement vectors are computed without pref- 
erence in azimuthal direction. We also multiplied the cross- 
correlation functions by a mask d(x\y') so that the maximum 
of the cross-correlation function is forced to be within a distance 
of Cs,\im = 12 pixels from its center. The typical time interval 
between consecutive images is in the range from 60-90 s, i.e., a 
feature moving at the photospheric sound speed of 8 km s"^ 
would travel 480-720 km corresponding to 6-9 pixels. This jus- 
tifies our choice of c^jim, which also takes into account some 
numerical errors. The position of the maximum of the cross- 
correlation function is calculated with subpixel accuracy by a 
parabolafit to the neighboring pixels. The numerical accuracy of 
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the parabola fit is about one fifth of a pixel or 16 km on the solar 
surface, which corresponds to about 200 m s"^ for proper mo- 
tions measured from a single pair of G-band images. Therefore, 
many flow maps have to be averaged to determine reliable hori- 
zontal proper motions. 

3.4. LCT data products 

Once the individual flow maps are calculated they are saved in 
binary format. In addition, average maps of horizontal speed and 
flow direction as well as the x- and j-components of the horizon- 
tal flow velocity (V;^, "^y) are stored in native IDL format. Some 
auxiliary variables are saved as well so that they can be used, 
e.g., in annotating plots depicting the flow fields. A sample of 
such plots for a 60-minute average flow field is shown in Fig. 2. 
In Fig. 2a horizontal proper motions are plotted as red arrows 
with a 60-minute averaged G-band image as a background. The 
moat flow starting at the sunspot penumbra and terminating at 
the surrounding G-band network is clearly discernible. 

In Fig. 2b we used an adaptive thresholding algorithm to dis- 
cern between granulation, G-band bright points, and strong mag- 
netic features. Indiscriminately, we used a fixed intensity thresh- 
old of /mag =0.8 for strong magnetic features and an adaptive 
threshold for G-band bright points, which can be given as 

/bp = 1.15 + 0.2(1 -yu) , (4) 

where = cos(^) is the cosine of the heliocentric angle 6. 
The darkest parts of sunspots (umbrae) and pores can be iden- 
tified using another fixed threshold of /dark = 0.6, while sunspot 
penumbrae cover intermediate intensities from /dark to /mag, al- 
lotting the range /mag to /bp to granulation. The adaptive thresh- 
old was necessary, as a first order approximation, to account for 
the center-to-limb variation (CLV) of the G-band bright points, 
which exhibit much higher contrasts near the solar limb. This 
adaptive thresholding algorithm allows us to study the proper- 
ties of horizontal proper motions for diflTerent solar features. 

Apart from the conventional way of displaying velocity vec- 
tors as arrows we present two-dimensional high-resolution speed 
(Fig. 2c) and azimuth (Fig. 2d) maps. In these maps, the physi- 
cal quantities are computed for each individual pixel so that the 
fine structure of the flow field becomes accessible. The color 
scale for the speed values is the same for all plots in this study. 
Indeed, we used the same color scale for all flow maps in the 
database so that flows for diflTerent scenes on the Sun can be di- 
rectly compared. In the azimuth map the direction is encoded in 
a 12-color compass-rose, which can be found to the very right of 
this panel. In principle more colors could be used to illustrate the 
flow direction. However, such plots would become very crowded 
and are very hard to interpret. The essential features of the flow 
field, i.e., inward motion of the penumbral grains in the inner 
penumbra and outward motions related to Evershed and moat 
flows, can easily be identified. 

Divergence (Fig. 2e) and vorticity maps were also compiled 
for each time-series. The maximum values for divergence (1 .5 x 
10"^ s"^) and vorticity (0.6 x 10"^ s"^) in these maps are an order 
of magnitude higher than values found for quiet Sun granulation 
in previous studies. However, higher values are not surprising, 
since we calculated them for each pixel in the high-resolution 
maps capturing more of the small-scale motions. Reassuringly, 
our values of the 10* percentile for divergence (4.3 x 10"^ s"^) 
and vorticity (1.9 x 10""^ s"^) are essentially the same as reported 
by Simon et al. (1994). On the other hand, in the present study 
the divergence is four times and the vorticity is two times higher 



than the values presented by Strous et al. (1996) but these values 
were computed in the proximity of an emerging active region. 
In summary and keeping in mind the diflTerent LCT parameters 
used in the diflTerent studies, the statistical properties of the flow 
fields as presented in this study are in agreement with previous 
investigations. 

In order to visualize the temporal evolution of the horizontal 
proper motions, we computed forward (Fig. If) and inverse cork 
maps. Evenly spread test particles are allowed to float forward in 
time with a given horizontal speed for a certain time interval (see 
Molowny-Horas, 1994). In the inverse cork map the particles 
float backward in time. This is accomplished by simply revers- 
ing the sign of the velocity components. The forward cork map 
is used to visualize regions of converging flows and the inverse 
map is a good tool to study divergence regions. We tracked test 
particles for consecutively two, four, and eight hours. These par- 
ticles were initially distributed on an equidistantly spaced grid 
with a spacing of 10 pixels, i.e., one particle was placed every 
0.8 Mm. The most conspicuous feature of the forward cork map 
are the tracer particles outlining the network of G-band bright 
points, which corresponds to the supergranular boundaries . 

We prepared overview web pages for the respective dates, 
when suitable time-series of G-band images were available, 
which contain all the six plots of Fig. 2 along with vorticity and 
inverse cork maps. The results of this study will ultimately be 
published as a small GAVO^ project as a value-added product of 
the Hinode database. 

3.5. Timing issues related to the image capture 

LCT delivers localized displacements observed in image pairs. 
These displacements in conjunction with the time, which has 
elapsed between the acquisition of both images, yield localized 
velocity vectors. Therefore, accurate knowledge of the time in- 
terval between consecutive images used in the LCT algorithm 
is essential. The time interval for the high-cadence image se- 
quence observed on 2007 June 4 has a bimodal distribution with 
values ^t of 14.4 s (60.7%) and 16 s (39.3%). The average value 
is A/" ^ 15.0 s. The diflTerence of 1.6 s is an artifact of the polar- 
ization modulation. The polarization modulation unit (PMU) is 
located just behind the telescope exit slit within the optical tele- 
scope assembly but before the tip-tilt mirror employed by the 
correlation tracker (see Tsuneta et al., 2008). A common CCD 
camera is assigned to both the broad- and narrowband filter im- 
agers (BFI and NFI). The critical timing between camera and 
PMU is handled by the focal plane package. The PMU is a con- 
tinuously rotating waveplate, which is always turned on - even 
for non-polarimetric data such as G-band images. Its rotation pe- 
riod is r = 1.6 s and all exposure timing is controlled with the 
clock of the PMU. This is the reason for the non-uniformity in 
the observed time-interval A^. 

However, we do not find a bimodal pattern in the LCT dis- 
placements but only fluctuations related to evolving features 
on the Sun and some residual numerical eflfects. Therefore, we 
opted to use the fixed time interval M in the data analysis. 
Sometimes a 'traflftc jam' in the data transfer might result in even 
largertiming errors. On the other hand, averaging individual LCT 
maps over an hour (or longer) will only result in velocity errors 
of less than a tenth of a precent, i.e., the speed measured by the 
LCT algorithm is not significantly aflfected. In summary, accu- 
rate timing has to be ensured to obtain reliable LCT flow maps 
and all data were checked for consistency between recorded time 

^ www.g-vo.org 
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Fig. 2. (a) The average (60-minute) G-band image after correction of geometrical foreshortening corresponding to the G-band 
image shown in Fig. la. The red arrows indicate magnitude and direction of the horizontal proper motions. Arrows with a length 
corresponding to the grid size indicate velocities of 1 km s"^. (b) Adaptive and fixed intensity thresholds are used to identify 
conglomerates of G-band bright points, granulation, penumbrae, umbrae, and pores. Color codes are used to point out (c) flow 
speed and (d) direction in high-resolution flow maps. The flow direction is encoded according to the 12 colors of the compass rose. 
Sources and sinks in a flow field can be identified in (e) a divergence map. The (f) forward cork maps provides further means of 
visualizing converging motions. Test particles are displayed after they have followed the flows for 2 (blue), 4 (orange), and 8 h (red). 
The conspicuous network of corks is related to the spatial scales of the meso-and supergranulation. 



stamps and measured horizontal displacements. Since the indi- 
vidual flow maps during the first and last two minutes of the 
one-hour sequences do not reflect the true proper motions but 
are artifacts of the subsonic filtering, we excluded them from the 
calculation of the average flow maps. 



4. Results 

4.1. Statistical properties of flow maps and time cadence 
selection 

The one-hour time-series on 2007 June 4 contains 238 G-band 
images, i.e., the time cadence is 15 s. This higher temporal 



resolution allows us to study the intrinsic accuracy of the LCT 
algorithm. We calculated LCT maps using seven diflferent time 
intervals A^ = 15, 30, 60, 90, 120, 240, and 480 s. Note that the 
time interval over which flow maps are averaged is reduced to 
Ar = 3600 s - A^, i.e., in case of the longest time cadence by as 
much as 8 min. However, these slightly diflferent averaging times 
will not change the results discussed below. When using all in- 
dividual LCT maps to compute the average horizontal proper 
motion we refer to such data as the entire sequence. On the other 
hand, when we split the entire sequence into four disjunct sets, 
we refer to them as interleaved data sets, i.e., every fourth LCT 
map is employed to compute the average horizontal proper mo- 
tion. Since these flow maps cover exactly the same period of 
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Fig. 3. Mean horizontal flow speed as a function of averaged 
LCT maps for time cadences of 15 s (thick dashed), 30 s 
(thick dash-dotted curve), 60 s (thick dash-triple-dotted curve), 
90 s (thick long dashed curve), 120 s (thin dashed), 240 s (thin 
dotted), and 480 s (thin soHd). 



time, diflTerences can be directly attributed to the numerical ac- 
curacy of the LCT algorithm. 

For the entire sequence and interleaved data sets we com- 
puted statistical parameters describing the distribution of hori- 
zontal proper motions for granulation in the vicinity of active 
region NOAA 10960. We used the adaptive thresholding al- 
gorithm (Eqn. 4) to select only granulation excluding G-band 
bright points. The proper motions, thus, refer to plasma motions 
in the absence of any strong magnetic field concentrations. Such 
selection facilitates comparing horizontal flow speeds and their 
distributions in all the cases of the present work. 

We calculated the mean v, median Vmed, maximum Vmax, and 
10* percentile vio of the horizontal flow speeds. Since the maxi- 
mum speed Vmax is only based on a single value, it is easily influ- 
enced by numerical errors and the data calibration. The 10* per- 
centile vio is more robust because it describes a property of the 
entire distribution, i.e., the high velocity tail. Along with these 
quantities we also calculated the variance cr^, standard deviation 
(Ty, skewness yi^y, and kurtosis 72,v The last two statistical pa- 
rameters describe the deviation of the distribution from a normal 
distribution. 

We find an increase of the average velocity with increasing 
time cadence Ar starting from about 0.40 km s"^ for Ar = 15 s, 
arriving at a maximum value of about 0.47 km s"^ for A^ = 60- 
90 s, and then decreasing from about 0.45 km s"^ for A^ = 120 s 
to 0.34 km s"^ for Ar = 240 s reaching the lowest value of about 
0.23 km s"^ for Ar = 480 s. Other statistical parameters such as 
Vmed, Vio, Vmax, CTy' foUow the samc trend. 

The initially increasing values can be explained by the time 
required for a solar feature to move from one to the next pixel. 
Three velocity values have to be considered: (1) the photo- 
spheric sound speed of ^^8 km s"\ (2) the maximum photo- 
spheric velocity of '^2 km s"^ measured by LCT techniques, and 
(3) the average speed for the proper motion of the granulation of 
^^0.5 km s"^. For A^ = 15 s the average displacement is around 
one tenth of a pixel, whereas the numerical accuracy for a single 
measurement is only one fifth of a pixel, because the maximum 
of the cross-correlation function can only be determined with 



this precision. Thus, in this case a solar feature has not suflftcient 
time to move, which results in underestimating its velocity. 

In case of A^ = 60-90 s, the horizontal displacement is suf- 
ficiently large so that a feature could have moved to one of the 
neighboring pixels. The speed in individual LCT maps is now 
well within the range, where numerical accuracy issues are neg- 
ligible. Starting at Ar = 120 s the mean velocity becomes lower, 
while there is suflftcient time for a solar feature to move quite 
some distance, the feature might have evolved too much, so that 
the LCT algorithm might not any longer trace the same feature. 
This leads to diminished horizontal velocities. 

Thus, 60-90 s is the good choice for measuring of horizontal 
flow speeds with LCT techniques. In this range of the time ca- 
dence, the mean values v of the interleaved data sets are essen- 
tially the same. Their deviations are much smaller than the pre- 
viously discussed systematic trends. In summary, all flow maps 
for the database were calculated using Ar = 60-90 s. If the time 
interval A^ was shorter, a multiple of the time interval Ar' = n/S.t 
with n = 2, 3, or 4 was used. 

4.2. Determining the duration of the time averages 

How many individual LCT maps have to be averaged to yield 
a reliable flow map? As previously discussed, the parabola fit 
to the maximum of the cross-correlation sets one limitation. 
However, there are also method-independent issues to be con- 
sidered. Solar features evolve over time so that a global pattern 
reveals itself only after averaging many individual LCT maps. 
We computed the mean horizontal flow speed as a function of 
the number of individual LCT maps that were used to arrive at 
an average flow map. The number of flow maps corresponds to 
the time interval AT over which the individual flow maps were 
averaged. 

Figure 3 presents this functional dependence for the time 
cadences from A^ = 15 s to 480 s. All curves start with high 
velocities when only a few individual LCT maps are averaged. 
It takes about 20 min before the curves level out and approach 
an asymptotic value indicating that these average flow maps are 
still dominated by the motions of fine structure contained within 
the sampling window. As discusses in the context of the statis- 
tical properties of the flow maps, short time cadences A^ tend 
to underestimate the flow speed. If flow speeds have be com- 
puted for time intervals shorter than 20 min, in particular for 
small-scale features, feature tracking methods are more appro- 
priate than LCT techniques. 

We find that all curves up to A^ = 90 s are stacked on top of 
each other without crossing the next higher curve at any point. 
Starting at A^ = 120 s we find that the curves corresponding 
to the longer time cadences cross the other curves after about 
20-25 min. This is another indication that solar features have 
evolved too much so that LCT fails to properly track their mo- 
tion. This behavior provides an explanation for the spread of ve- 
locity values found in literature. In particular, short time- series, 
as often encountered in ground-based observations, might be bi- 
ased towards higher velocities. In summary, our choice of 60- 
minute averages for LCT maps is a conservative one giving solar 
features sufl^icient time to reveal the global flow pattern. 

4.3. Selection of the sampling window 

How do the horizontal proper motions depend on size and 
FWHM of the Gaussian kernel used in the LCT algorithm? To 
answer this question, we calculated horizontal proper motions 
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Fig. 4. Horizontal flow speed {left) and azimuth maps {right) 
measured using Gaussian kernels with 64 x 64 pixels and var- 
ious FWHM. The FWHM of 7.5, 15, 22.5, and 30 pixels cor- 
respond to 600, 1200, 1800, and 2400 km on the solar surface, 
respectively. The bottom row depicts smoothed versions of the 
FWHM = 1200 km speed and azimuth maps. The color coding 
is the same as in Fig. 2. 



using a Gaussian kernel with 64 x 64 pixels, which is equivalent 
to 5120 km X 5120 km on the solar surface. This larger ker- 
nel was chosen to encompass successively broader FWHM. We 
choose four FWHM of 7.5, 15, 22.5, and 30 pixels corresponding 
to 600, 1200, 1800, and 2400 km, respectively. The FWHM of 
1200 km matches the size of a granule. Individual LCT maps are 
produced from image pairs separated by 60 s in time. We com- 
puted the statistical parameters relating to granulation for the 
entire sequence and the interleaved data sets. All statistical pa- 



rameters are decreasing with increasing FWHM. For FWHM = 
600 km the mean velocity is v = 0.48 km s"\ which decreases 
to V = 0.47 km s"^ for FWHM = 1200 km, v = 0.41 km s"^ 
for FWHM = 1800 km, and further to v = 0.37 km s"^ for 
FWHM = 2400 km. Small-scale features, which exhibited the 
highest proper motions, are lumped together with regions of low 
flow speeds when the FWHM increases. This eflTect is displayed 
in Fig. 4, which shows the speed and azimuth maps for the four 
FWHM. As before, for a given FWHM all statistical parameters 
agree with each other for up to three significant digits, i.e., there 
is no indication that the algorithm's numerical accuracy depends 
on the FWHM of the kernel used in LCT. 

In the case of FWHM = 1200 km, the statistical parameters 
are very close to the ones calculated for the same FWHM but 
with a kernel of 32 x 32 pixels. This is not surprising, since the 
Gaussian kernel assigns a much stronger weight to features in the 
center of the subimage, i.e., the periphery in the 64 x 64 pixels 
FOV has only a small influence in determining the displacement 
vector for a pair of subimages. This of course only holds true as 
long as the wings of the Gaussian do not significantly extend be- 
yond the edges of the kernel. Convolving the average flow map 
(FWHM = 1200 km) with a Gausssian kernel, which had a size 
of 32 X 32 pixels and FWHM = 26.4 pixels (2112 km), we ar- 
rive at a smoothed version (see bottom row in Fig. 4), which is 
virtually identical to the flow map with FWHM = 2400 km. 

In conclusion, it makes no diflTerence, if one uses a larger 
FWHM while computing LCT maps or one smoothes the maps 
after computation. In both cases, the results are virtually the 
same. Since only minor changes in the LCT results were ob- 
served for kernels with 32 x 32 pixels as compared to 64 x 64 
pixels, the smaller kernel was chosen significantly reducing the 
computing time. Furthermore, the smallest FWHM produces 
the most detailed flow maps. However, we chose a FWHM of 
1200 km favoring the spatial scales of granulation, which covers 
the largest fraction of the observed area. Additional smoothing 
can still be applied in the later data analysis stages to either re- 
duce noise or to track flows on larger spatial scales. For case 
studies of sunspots fine structure smaller FWHM might be more 
appropriate. 

4.4. Numerical errors in calculating flow maps 

We computed the pixel-to-pixel rms-error for the magnitude and 
direction of the flow velocity using the interleaved data sets for 
diflferent time cadences and for diflTerent FWHM of the sam- 
pling window. Six diflference maps (sets 1-2, 1-3, 1-4, 2-3, 
2-4, and 3-4) can be computed from the four interleaved data 
sets, thus for each pixel we can derive the errors, which are pri- 
marily due to numerical errors inherent to the LCT algorithm . In 
case of diflferent time cadences, the rms-error in velocity is 15- 
90 m s"^ from shortest to longest cadence. The corresponding 
rms-error in direction is 5°-30°. However, for cadences in the 
range of 60-90 s, which is the range used to create the database, 
the rms-error of the velocity is typically in the range from 3 5- 
70 m s"\ while the values for the direction vary by as much as 
10°-15°. The largest variations in direction are observed near 
the boundaries of patches showing coherent flows. Even after 
correcting the In ambiguity in the diflference maps, we find high 
values at these locations. As a side note, the rms-error in di- 
rection justifies our choice of a color wheel with only twelve 
segments in the display of the azimuth maps. One segment cov- 
ers 30° so that pixel-to-pixel variations of about ±15° are sup- 
pressed. Otherwise, the azimuth maps would appear too crowded 
distracting from the overall flow pattern. 
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Fig. 5. Time-averaged G-band images with horizontal flow vectors computed for Ar = 1, 2, 4, 8, and 16 h. Arrows with a length 
corresponding to the grid spacing indicate velocities of 1 km s"^. The background images are displayed in the intensity range 
0.5 <///o < 1.5. 



In the case of diflTerent FWHM of the sampling window, rms- 
errors for speed and direction decrease with increasing FWHM. 
The velocity error is 35 m s"^ and the error in direction is about 
10° for a FWHM of 1200 km. These errors decreases to 15 m s"^ 
and 5° for a FWHM = 1800 km and further to 10 m s'^ and 4° 
for a FWHM = 2400 km. Here, the decreasing rms-errors can 
be attributed to the smoothing eflfect of a wider sampling win- 
dow, i.e., more pixels are used with higher weights in the cross- 
correlation. Furthermore, the rms-error in magnitude and direc- 
tion for the time cadence Ar = 60 s is nearly the same regardless 
of the size of the Gaussian kernel (32 x 32 pixels vs. 64 x 64 
pixels. 

Finally, we calculated Pearson's correlation coeflftcient be- 
tween LCT maps of the interleaved data sets. Pearson's correla- 
tion coeflftcient indicates the degree of a linear relationship be- 
tween two variables. A positive value of unity indicates that the 
data sets are identical disregarding a linear scaling factor. The 
linear correlation coefl^icient for diflTerent time cadences A^ de- 
creases from 0.99 to 0.93 starting at the shortest and ending at 
the longest cadence. The high degree of correlation indicates that 
all essential features of the flow field are capture by the LCT 
algorithm. The monotonic decrease, however, indicates that nu- 
merical errors increase when the cadences A^ become too large. 

4.5. Long-lived features In flow maps 

The one-hour time interval over which the LCT maps are aver- 
aged is not sufi&cient to identify features, which need longer to 
evolve such as meso- and supergranulation. Therefore, to clearly 
identify the boundaries of these large-scale convective cells and 
to visualize the eflTects of longer time averages, we averaged LCT 
maps over Ar = 1, 2, 4, 8, and 16 h utilizing the long-duration 
data set with 960 images and a time cadence A^ = 60 s. 

As before, we calculated the statistical parameters for the 
entire sequence and the interleaved data sets for granulation. 
All velocity values decrease with increasing time intervals over 
which, the individual LCT maps are averaged. For Ar = 1 h 
the mean velocity is v = 0.44 km s"^ which decreases to 

V = 0.38 km s"^ for AT = 2h,v = 0.34 km s"^ for AT = 4h, 

V = 0.30 km s"^ for AJ = 8 h, and further to v = 0.23 km s"^ for 
Ar = 16 h. The mean speed approaches a value for the global 



flow field with increasing AT. However, the value for AT = 1 h 
is only slightly smaller than previously computed for the high- 
cadence sequence. Such small deviations (< 0.05 km s"^) reflect 
only minute diflferences between the scenes on the solar surface. 

Fig. 5 contains the time-averaged G-band images with su- 
perposed arrows indicating speed and direction of the horizon- 
tal proper motions for Ar = 1, 2, 4, 8, and 16 h. The region 
shows granulation and G-band bright points to the west of ac- 
tive region NOAA 1 1930. Magnetic features dominate the long- 
duration time averages, i.e., the flow speed is low where strong 
magnetic fields are present in the chromospheric network and 
conversely the speed is high where G-band bright points outline 
the boundaries of large-scale convective cells. The high speeds 
associated with the local convective pattern of granules have 
diminished for the long-duration flow maps and only the con- 
verging motion towards the cell boundaries remains explaining 
the statistical properties of the velocity values discussed above. 
The overall visual impression of the vector maps is that the ar- 
rows are more ordered in the long-duration maps, whereas in the 
short-duration maps (AT = 1 and 2 h) a larger scatter of the 
flow vectors is observed on smaller scales. Nonetheless, the im- 
print of the the meso- and/or supergranulation is already visible 
in the short-duration flow maps and becomes more prominent 
the longer the time interval over which the LCT maps are aver- 
aged. Strong converging motions can be found in Fig. 5 near the 
vertical alignment of G-band bright points in the north-east cor- 
ner of the FOV and at the boundary of the larger supergranular 
cell in the south-east corner of the FOV. The supergranule also 
contains substructure on smaller scales, i.e., a strong divergence 
center exactly in the central FOV, which can be clearly seen after 
averaging for at least Ar = 2 h. 

We plotted relative frequency distributions of flow fields 
spanning time intervals of Ar = 1, 2, 4, 8, and 16 h in Fig. 6 
to gain insight into the statistical properties of the long- duration 
data sets - this time both for granulation and for G-band bright 
points. The statistical parameters characterizing these distribu- 
tions are provided in Tab. 1 for reference. The distribution de- 
scribing granulation for the shortest time interval (Ar = 1 h) is 
the broadest and has an extended high velocity tail. Interestingly, 
for velocities up to 0.6 km s"^ this distribution is virtually the 
same as the distribution for an averaging time, which is twice 
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Table 1. Parameters describing the horizontal proper motions of granulation and G-band bright points calculated over Ar 
1, 2, 4, 8, and 16 h. 
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Fig. 6. Relative frequency distributions of horizontal proper motions for quiet Sun {left) and G-band bright points {right) calculated 
over Ar = 1 h (solid curve), 2 h (dotted curve), 4 h (dashed curve), 8 h (dash-dot-dotted curve), and 16 h (dash-triple dotted curve) 



as long (Ar = 2 h). The only difference is the high velocity 
tail. This indicates that proper motions on small scales still make 
their presence known, if individual LCT maps are not averaged 
for at least two hours. The longer Ar, the more are the peaks of 
these distributions shifted towards lower velocity values (from 
0.43 to 0.23 km s"^) and their standard deviations are progres- 
sively becoming smaller (from 0.24 to 0.12 km s"^). The pro- 
gression of the frequency distributions shown in the left panel 
of Fig. 6 supports the conclusion that the essential features of 
long-duration LCT maps have been captured for Ar = 8-16 h. 

The frequency distributions for G-band bright point differ in 
some aspects from those for granulation. The high velocity tail is 
less prominent and all statistical parameters describing the dis- 
tributions are reduced by about a factor of two. The two distribu- 
tion with the shortest time intervals (Ar = 1 and 2 h) show a hint 
of a bimodal distribution, and they are skewed towards higher 
velocity values. However, this might be an artifact of the adap- 
tive thresholding algorithm, since the areas covered by G-band 
bright points are smaller and well defined in the short-duration 
intensity maps. Thus, considering the FWHM = 1200 km of the 
sampling window, a larger contribution from granulation is ex- 
pected, if more isolated G-band bright points are present in the 
maps, which are used for threshholding. It should be noted that 
Eqn. 4 was slightly modified to accommodate the longer time in- 
tervals, which leads to a fuzzier appearance of the area covered 
by G-band bright points and results in a diminished contrast of 
the G-band bright points. Since no contemporary magnetograms 



with comparable spatial resolution were available, we cannot 
comment on the influence of flux emergence or dispersal during 
the observed time interval. However, active region NOAA 10930 
showed pronounced activity. In particular, the penumbra of the 
small sunspot just to the east of the FOV decayed and resulted in 
continuous flaring in the active region. 

Even though long-duration time averages are an important 
tool when studying large-scale convective patterns or the per- 
sistent motions in an active region, the scarcity of such data 
sets argues against their use for comprehensive and compara- 
tive studies. Since most of the characteristics of flow fields are 
already captured in one-hour averages, we opted for Ar = 1 h. 
Furthermore, selecting Ar = 1 h allows us to study changes 
in long-duration time- series by computing averaged flow maps 
every 30 min. Another consideration is that one hour is more 
then ten times the typical lifetime of granules so that the proper 
motions of individual granules should be negligible and global 
motion patterns will reveal themselves. 

4.6. Frequency distributions for different solar features 

The simplest approach to describe flow fields would be to com- 
pute the overall frequency distributions for a particular FOV. 
However, this simplistic approach is insufficient to recover the 
underlying physics of plasma motions in the presence (or ab- 
sence) of strong magnetic field. We used the adaptive threshold- 
ing algorithm (Eqn. 4) described in Sect. 3.4 to compute fre- 



10 M. Vema and C. Denker: Horizontal flow fields observed in Hinode G-band images 

Table 2. Statistical parameters describing the frequency distributions of the horizontal proper motions of various solar features. 
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quency distributions for granulation, penumbrae, umbrae/pores, 
G-band bright points, and the entire FOV regardless of the fea- 
tures contained in this region. For this case study, we used the 
high-cadence sequence (see Fig. 2b for the thresholded image) 
with = 60 s, Ar = 1 h, and a FWHM = 1200 km. The re- 
spective distributions are shown in Fig. 7 and the corresponding 
statistical parameters are summarized in Tab. 2. Similar plots and 
tables will be included in the database subsuming the more than 
200 data sets, which have been analyzed as part of this study. We 
provide Fig. 7 and Tab. 2 to facilitate comparison with other case 
studies. However, these data are not representative (in the sense 
of a mean value) for all data sets contained in the database. 

A barely detectable shoulder in the frequency distribution for 
the entire FOV and the extended high velocity tail hint already 
that this distribution contains contributions from various solar 
features. Its mean velocity v = 0.44 km s"^ is slightly smaller 
than the corresponding value for granulation v = 0.47 km s"^ 
which dominates the FOV. The distribution for granulation is 
broader and a low value of kurtosis ij2,v = 0.08) leads to a flat- 
ter peak, where any indication of a shoulder is absent. There is 
a noticeable diflTerence in the distributions of strong magnetic 
elements and granulation. The distributions for G-band bright 
points, umbral and penumbral regions are narrow, have sharp 
peaks, and are shifted towards lower velocities. The mean veloc- 
ity for these regions varies from v = 0.30 km s"^ for penum- 
brae to V = 0.23 km s"^ for umbrae/pores and G-band bright 
points. Interestingly, the distributions for umbrae/pores and G- 
band bright point are virtually identical, while that for penum- 
bra has significant contributions at velocities above 0.4 km s"^. 
This is indicative of the more complex flow fields in penumbrae, 
where penumbral grains move preferentially in the radial direc- 
tion - inward in the inner penumbra and outward in the outer 
penumbra. This also illustrates that some of the small-scale hor- 
izontal proper motions can be capture with the current imple- 
mentation of the LCT algorithm. 



4. 7. Flow maps for different spatial resolution 

Finally, we would like to address the question on how the spatial 
resolution aflTects the determination of the horizontal proper mo- 
tions. We used the high-spatial resolution sequence of Sect. 2.3 
and treated it exactly the same as all the other data with the ex- 
ception that the G-band images were sampled at 40, 80, 120, 
160, and 200 km after correction of geometrical foreshortening. 
Multiple of 40 km were chosen to match the Hinode SOT/BFI 
pixel size of 0.055''. Obviously the number of pixels in the sam- 
pling window had to be adjusted. However, shape, size, and 
FWHM = 1200 km of the sampling window were not changed. 




Velocity [kms-^] 

Fig. 7. Relative frequency distributions of horizontal proper mo- 
tions computed for all solar features (solid curve), G-band bright 
points (dotted curve), granulation (dashed curve), umbra (dash- 
dotted curve) and penumbra (dash-dot-dotted curve). 



The values describing the respective frequency distributions 
are given in Tab. 3. They were computed for areas with granu- 
lation covering the full FOV However, to visualize the minute 
changes in the flow maps, we show in Fig. 8 only an area of 
8 Mm X 8 Mm. In the top row of Fig. 8 flow vectors are super- 
posed on one-hour average G-band images with diflTerent image 
scales. This scene on the Sun is dominated by converging mo- 
tions towards the X- shaped alignment of G-band bright points. 
The grid spacing corresponds to 320 km and the length of the 
arrows was chosen so that an arrow with a length of exactly 
the grid spacing corresponds to 0.5 km s"^ The diflTerence are 
so minute that they only show up in diflTerence images of the 
LCT maps. The bottom row of Fig. 8 shows the flow speed for 
each pixel in the FOV at the same color scale as used in all 
the other figures. The overall appearance to the flow field is the 
same. However, the low resolution maps look blockier due to the 
coarser sampling and some of the fine structure starts to fade out. 

In summary, the average velocity diminishes from 
0.54 km s"^ at the highest spatial resolution to 0.47 km s"^ at the 
lowest resolution. This trend is the same for all other parameter 
with the exception of kurtosis and skewness, which show some 
(negligible) scatter. Changes of less than 15% in velocity cannot 
explain the broad range of velocity values for horizontal proper 
motion reported in literature. Note, however, that Hinode data is 
not susceptible to the adverse aflTect of seeing, i.e., ground-based 
LCT measurements will be much more affected depending on 
the spatial resolution. Even though seeing should not introduce 
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Table 3. Statistical parameters describing the frequency distributions of the horizontal proper motions for various image scales. 
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Fig. 8. Time averaged G-band images with horizontal flow vectors {top) and flow speed maps (bottom) for diflferent spatial resolu- 
tion. Arrows with the length corresponding to the grid spacing indicate velocities of 0.5 km s"^ 



a systematic bias in LCT (see November & Simon, 1988), it feature, location on the Sun, or solar activity. In the following 
will still aflTect the noise in the LCT measurements. list, we summarize our choice of LCT parameters. 



5. Conclusions 

Many case studies exist in the literature describing horizontal 
proper motions based on LCT or FT techniques. Even though 
most of them agree on the morphology of the observed flows, 
significant diflTerences are found when quantifying the flow prop- 
erties. Besides obvious diff'erences inherent to the techniques, 
the choice of parameters such as sampling window, time ca- 
dence, and duration can significantly impact the outcome. Some 
results of previous studies are provided in Tab. 4 for convenience 
and to ease comparison with the present investigation. 

We have presented the implementation of an LCT algorithm, 
which was used to create a database of flow maps derived from 
time-series of G-band images observed with Hinode/SOT. The 
parameter study and error analysis will also be beneficial to other 
studies using LCT techniques. Even for observations from the 
ground our results provide guidance, since LCT techniques are 
not biased by seeing (see November & Simon, 1988) so that our 
error estimates can be understood as a lower limit. 

Justifying the choice of parameters for LCT and FT algo- 
rithms is always a challenging task, which should be driven by 
the scientific purpose of the study. In the present study, the em- 
phasis was on creating a database of flow maps, which can be 
used in statistical investigations regardless of the type of solar 



□ The flow maps are based on time-series of G-band im- 
ages with cadences between 60 s and 90 s. If the ca- 
dence is shorter, features with low velocities cannot be accu- 
rately tracked, whereas in longer cadences features will have 
evolved too much so that the algorithm might not recognize 
them any longer. Our cadence selection is conservative in the 
sense that we limit our database in favor of better compara- 
bility. Note that there is a significant number of G-band time- 
series with cadences of about 2 min, which are not included 
in our database. 

□ The evolution of individual features (granules, bright points, 
penumbral grains, umbral dots, etc.) dominate flow fields on 
short time scales. Therefore, averaging over time scales sig- 
nificantly longer than the lifetime of the aforementioned fea- 
tures is necessary to yield the global flow field. Our choice of 
Ar = 1 h over which the flow maps are averaged ensures that 
the global flow fields has emerged from the motions of indi- 
vidual features and that sufl&cient flow maps were averaged 
to reduce the numerical rms-errors for magnitude and direc- 
tion of the flow vectors to reasonable values of 35-70 m s"^ 
and 10-15°, respectively. Note that Ar = 1 h is not an appro- 
priate choice for studies focusing on meso- and supergranu- 
lation because the associated flow pattern is still very noisy 
and could be more easily perceived in longer time averages. 
However, long-duration time- series are rare to facilitate such 
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Table 4. Summary of LCT results in previous studies. 



a FWHM 


dgrid 




AT 




Vmed 


^max 


Remarks 


[pixel-i] 




[s] 


[min] 


[km s-i] 


[km s-i] 


[km s-i] 




November & Simon (1988) and November (1989) 

0.250'' 33" 2.0" 67 80 0.5-1.0 (for source regions) 


continuum images at /1517.5 nm, Universal Birefringent 
Filter (UBF), 12 s exposure time, Dunn Solar Tele- 
scope/Sacramento Peak 


Brandt et al. (1988) 
0.035" 2.4" 


0.8" 


60 


79 


0.67 




1.2 


broad-band (5.4 nm) images at /1469.6 nm, Solar 
Optical Universal Polarimeter (SOUP), 20 ms exposure 
time, Swedish Solar Vacuum Telescope/Observatorio 

UCl XvUL|UC UC lUb iVlULllaLllUJ> 


Title et al. (1989) 
0.161" 4.0" 


2.5" 


60 


28 


0.37+0.45 




<1.2 


broad-band (100 nm) images at A600 nm, SOUP, 
short-exposure images recorded on photographic film, 

vSnarplab 9 r>rnvirlp<s al<sn nronpr motion mpa<siirpmpnt<s 

for other FWHM down to 1.0" 


Bergeret al. (1998) 
0.083" 0.83" 


0.4" 


24 


70 


1.10+1.30 


0.70 


-4.0 


G-band images /1430.5 + 0.6 nm, SOUP, 20 ms exposure 
time, SVST, LCT parameter study and comparison of 
proper motions between granulation and network 


Shine et al. (2000) 
0.600" 4.8" 


2.4" 


60 


60 


0.49 


0.47 


1.5 


continuum images near the Nil /1676.8 nm line, 
Michelson Doppler Imager (MDI), Solar Heliospheric 
Observatory (SoHO), long-duration sequence of 45.5 h 



Note: a: image scale, FWHM: width of the sampling window, dgnd- grid spacing on which LCT flow vectors are computed. At: image 
cadence, AT: time interval over which individual LCT maps are averaged, v + cr^: average flow speed and standard deviation, Vmed- 
median value of the frequency distribution, and Vmax' highest observed flow speed. 



studies. Whenever, time- series with longer durations were 
available, we computed one-hour flow maps with an overlap 
of 30 min so that the temporal evolution of the flow field can 
be monitored, which is of particular interest for the investi- 
gation of explosive events such as flares, filament eruptions, 
and coronal mass ejections. 
□ In principle, the spatial resolution of Hinode/SOT would al- 
low to track features, which are smaller than one second of 
arc. Our choice of a Gaussian sampling window with 32 x 
32 pixels and a FWHM of 1200 km was again motivated by 
establishing a database of flow maps for statistical studies. 
Therefore, we used a FWHM, which corresponds approxi- 
mately to the size of a granule, which is one of the 'largest' 
elements of solar fine structures. Tracking flows on larger 
spatial scales can still be accomplished by smoothing the 
flow maps after the fact. 

In forthcoming studies, we will use the database of flow 
maps to study the statistical properties of pores, the motions in 
sunspot penumbrae, and their relation to the flow pattern ob- 
served in the moat of sunspots. Several years of G-band time- 
series and more than 1000 individual flow maps facilitate the 
study of such flows during the life cycle of solar features and en- 
vironment, i.e., as a function of solar activity or the complexity 
of the surrounding magnetic field. Once thoroughly tested, the 
value-added Hinode/SOT data will be made available as a small 
project within the scope of GAVO. 
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